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We experimentally verify uncertainty relations for mixed states in the tomographic representation 
by measuring the radiation field tomograms, i.e. homodyne distributions. Thermal states of single- 
mode radiation field are discussed in details as paradigm of mixed quantum state. By considering the 
connection between generalised uncertainty relations and optical tomograms is seen that the purity 
\ of the states can be retrieved by statistical analysis of the homodyne data. The purity parameter 

assumes a relevant role in quantum information where the effective fidelities of protocols depend 
critically on the purity of the information carrier states. In this contest the homodyne detector 
becomes an easy to handle purity-meter for the state on-line with a running quantum information 
protocol. 
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p ,■ I. INTRODUCTION 

' Quantum pure states can be faithfully described by their wavefunction ^ ix) . Experimentally producing and 
CJ^. measuring pure quantum states is impossible due to the imperfections both of the generation procedures and the 
measurement apparata. In particular quantum states of the optical fields are always mixed due to impossibility of 

(N ■ 

^ ' transmitting and detecting optical fields with 100% efficiencies. Density matrix p formalism encopasses the possibility 
0^ ■ of describing both pure and mixed state. In this context an important parameter is represented by the purity tt given 

<N ■ 

' by the trace of the squared density matrix operator with tt = 1 for pure states while tt < 1 for mixed ones. Moreover, 
TT assumes a relevant role in quantum information protocols where the fidelity, i.e. the rate of success, depends 
critically on the purity of the states, tt also can give a quantitative measure of the decoherence the pure state have 
suffered. 



Pure states satisfy the Schrodinger 



2j and Robertson [3| basic inequalities that generalize Heisenberg [J| uncertainty 
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• relations including contributions from covariance of conjugate quantum observables. In Ref. [5| (see also [6|) a new 
bound, higher than the Schrodinger- Robertson, which accounts for the contribution of the purity of mixed states was 
found. In Ref. Q the Schrodinger-Robertson uncertainty relation was expressed in terms of homodyne tomograms. 
The generic tomographic approach to quantum systems was reviewed in [8| . Combining the purity dependent bound 
and its expression in terms of homodyne tomograms give rise to a method for a simple determination of the state 



purity via homodyne detection developed in 9l-ll4 1 . 
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The problem of measuring the position and momentum was discussed also in connection with the tomographic 
approach recently in [l^. 

So far, evaluation of purity in continuous variable (CV) systems is obtained once the full state density matrix have 
been reconstructed by quantum tomographic methods. In this paper we prove that the purity of the analyzed state 
can be retrieved by a very simple and fast analysis of homodyne data thus giving the possibility of having an on-line 
monitor on the state. 

The aim of the present paper is twofold. On the one hand we express the purity tt of a mixed state in terms of 
homodyne tomograms, i.e. quantities amenable to an experimental determination. On the other hand, relying on the 
uncertainty relations of Ref. we derive a simple estimator for the purity of a thermal state which is independent 
from its to mog raphic expression. 



In Ref. [16| a probability representation of quantum mechanics was suggested in which states are described by 
standard probability distributions, called tomograms. This representation is based on the representation of the 
Wigner function W ( p, q) of a quantum state by means of the Radon integral transform or marginal distributions 



(optical tomograms) 



17 



I . These data representing the output of optical homodyne detectors allows reconstructing 
the Wigner function using the experimental data. This tomography procedure is nowadays a routine method for 
measuring quantum states (see for a review [l^. fioj). In the present paper, using ideas discussed in we apply the 
optical tomography to check generalized Schrodinger-Robertson uncertainty relations for conjugate quadratures in 
the case of a mixed quantum state, in particular, for measuring the effective temperature of a state. The idea of our 



consideration is based on suggestion [16[ (see, also [8|) that the homodyne tomogram is a primary object identified 
with the quantum state. Due to this, one can extract all information on the state properties including the purity from 
the measured tomogram only, avoiding the reconstruction of the Wigner function procedure. 
The paper is organized as follows. 

In section 2, starting from the symplectic forms we give an expression of the purity parameter and of the mean 
photon number of any mixed photon state in terms of measurable optical tomograms. For thermal states this 
expression provides the temperature of the field and the mean photon number thus suggesting to check the accuracy 
of homodyne detection by comparing photon statistics obtained in this way, i.e., from optical tomograms, and in 
independent photon counting experiment. 

In section 3, we review mixed state uncertainty relations to get an estimation of the purity independent of its 
tomographic expression, by saturating the bound set by the uncertainty relation found in [5|. Hence, measuring 
optical tomograms of photon states, we can evaluate the purity of the state, and, in the case of thermal states, study 
the dependence of the quantum bound on the field temperature, thus obtaining a sort of a thermometer for evaluating 
the temperature. Also, we obtain an estimation for the mean photon number. Finally, in section 4, an experimental 
comparison, for thermal states, is performed between our estimation and an independent measure of purity of the 
state, which allows for evaluating the accuracy of our approximate estimation. 
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II. TOMOGRAM AS PURITY AND THERMAL METERS 



A quantum state, either mixed or pure, is described by a density operator p with purity parameter tt given by: 

= Tr < 1 . 

According to {igI the state is described by a symplectic tomogram W (X, /i, , where X, /i, are real parameters, 
obtained by the Radon transform of the Wigner function W {p,q) (hereafter, h — 1): 



dpdq 
2n 



Tr 



p 6 [X - nQ - i^P 



(1) 



with S (^X — pQ — vP^ standing for 



X ^ pQ + vP 



Accordingly the generalized quadrature operator X (/i, v) depends parametrically on /i, v and its moments are given 
by: 



(Jl" (Ai, u)^ = I X"W {X, p,v)dX , n = 1, 2, 3, . 



like means and variances, in terms of homodyne quadrature statistics. In particular, 



X^ {p, v) = p^Q^ + v^P^ + 2pv 



( QP + PQ 



so that 



^ [x [p, v) - (x {p, v)y^ ^ = (Txx (m, i") = f^^<^QQ + f'^crpp + 2pvapQ 



(2) 



(3) 



Accordingly, the mixed state is characterized by the quadrature dispersions crQQ,app,apQ, which in turn can be 
determined by measuring axx {p, i^) for particular values of p, v. In case of optical tomograms, one has p = cos 9, v = 
sin 6* so that the optical field state is characterized by the field quadratures relative to 9 ~ 0,7r/2 and 9 = tt/4 : 



axx {^,0) , o-pp ^ crxx (0,1) , 

axx (^'^) -^Kx(l,0) + axx(0,l) 



O'PQ 

Analogously the photon number operator h — a^a — ^ i^Q^ + — 1 j is related to the moments 

{n) = ^[{x'il,0)) + {x'{0,l))-l_ . 



(4) 



Hence the accuracy of the homodyne detection could be assessed by comparing (h) obtained via optical tomograms 
with that measured by standard photon counting experiments. 

The purity tt of /5 is a functional of W {X, p, v) (see, for example, [8j) 

^ = Tr [p^] = 2^ / dXdYdpdv e'^^+^^yv (X, p, u) W (F, --p, ~v) 
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or equivalently 

I*2tt poo 

I I dXdY i 

2tt 



jT= 1 1 dXdY I de I dk 



^0 

having replaced W {X, fi, v) with the homodyne marginal distribution 

Wo {X, 9) = kW {kX, k cos 9, k sin 9) 

which is accessed in homodyne measurements. 

For Gaussian photon states W {X, ^, v) reduces to: 

/X, v) = ^ = exp 

y'2TTaxx (m, I') 

which inserted © into (O yields the well-known expression: 

1 



'^\I<^QQ(^PP - <^QP 



(6) 



(7) 



For a thermal state (cqq = <ypp = \ coth (^) ; ctqp = O) tt reduces to 



* = tanh ( ^ ) (8) 



with T measured in Kelvin x Kb/Tt-oj. 



III. MIXED STATE UNCERTAINTY RELATION 



The general uncertainty relation for a mixed state reads 



with TT the purity of the state p. The real, continuous and differentiable function $(7r), such that (f>(7r) > 1 in the 
interval < -ir < 1 , has the following piecewise analytic expression (extrema of intervals are given by 2(2fc+l)/3fc(fc+l), 
fc = l,2,...): 
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$(7r) = 2 - V27r - 1 - < tt < 1 



Besides, the function $(7r) can be approximated in the whole interval (0, 1) within 1% by the interpolating function 



$W = . (11) 
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FIG. 1: Relative difference between the bound and its approximation $^ given in Eq. $^ approximates $^ witliin a 

few percent s. 



In Fig. (IlJ we plot the relative difference between ^'^{tt), the square of the bound function, and <I>^(7r) in order to 
visualize how good is the approximation. 

We generalize the inequality (j9l) to arbitrary local oscillator phases 0, by using the tomographic uncertainty function 
F{6) introduced in Q|: 

1 



Then, by using $(7r) instead of $(7r), the uncertainty relation reads: 

1 



F{0)-- =Fi9) 



2V16 + - 18n^ 
81^2 



> , 



(12) 



(13) 



in this way, a relation between the tomographic uncertainty function F (9) and the state purity is set. 
We recall that F (9) is defined, by means of the variance of the homodync quadrature X 

r ~i 2 

axx{0)^ J X^WoiX, 9)dX - J XWoiX, 9)dX 

as: 



(14) 



c^xx 



(^+94(.xx(^)+.x.(^+f)) 



1 

4' 



We note that for = one has 



F (^)le=o '^QQ'^PP - ^QP 



(15) 



(16) 



So that, F {9)\g^Q > is exactly the Schrodingcr-Robertson inequality. Moreover, comparing the last expression with 
the thermal state purity given in Eq. ([7]) it is easy to see that: 

1 

The above expression of the uncertainty function F (9) in terms of tomograms was given in . 

The physical meaning of the function F (9) is the following. For a local oscillator phase 6* = 0, it is the determinant 
of the quadrature dispersion matrix, shifted by —1/4 as shown in Eq. (jl6|) . For a nonzero local oscillator phase 9, 



the function F (6) corresponds to the determinant of the dispersion matrix of the quadratures, which are measured 
in a rotated reference frame in the quadrature phase space. The non-negativity of the function F (9) imphes the 
fulfilhng of the Schrodinger-Robertson uncertainty relation for all the unitarily equivalent position and momentum 
operators, since the unitary transformations do not change the canonical commutation relations. The formula (|15|) 
simply expresses the determinant of the dispersion matrix for unitarily rotated position and momentum, in terms of 
tomographic probability distribution Wo {X, 9) . 

Thus, the tomogram must satisfy the inequality (jl2p . or (|13p . where the parameter tt is expressed in tomographic 
terms by Eq. ([5]). However, we can get an estimation of the purity in terms of the tomographic uncertainty function 
F (9) by saturating the inequality (|13p . In other words, we consider the minimum value F of the uncertainty function 
F {9) and estimate that for such a value the inequality is pretty near saturated. Then, by solving with respect to tt, 
we are able to express the purity as a function of F. This is a simple expression: 



2Vl + 4f 
2 + 9F 



(17) 



whose plot is shown in Fig. [2j 




FIG. 2: Purity vs. F the minimum value of F (9). F saturates the inequality ([13] 



As one can see. _F = corresponds to a pure state, tt = 1, and the purity is a smooth decreasing function, going to 
zero for F at infinity. 

On the other hand, for thermal states the purity parameter is directly related to the temperature, i.e. tt ~ 
tanh(l/2T), so that the previous inequalities and estimations can be translated in terms of temperature. So, in 
particular, for a thermal state the bound <f> in inequality (|13l) can be written in terms of the temperature T. 

Moreover, F (9) provides the temperature T {F) of a thermal state with the corresponding purity tt (F). We get 



T= [2tanh~^^] ^ 



2 tanh 



_i /2V1 + 4F 



(18) 



V 2 + 9F 

We recall that for a thermal state the tomographic uncertainty function F{9) does not depend on the local oscillator 
phase 9, as shown in so that F{9) = F, V^. 

Finally, we remark that for a thermal state the mean value of the photon number can be expressed as a function of 



the temperature as 



so we have a relation between our estimation of the temperature T (F) and the photon statistics. 

In other words, we may check the accuracy of our estimation of the purity and the temperature, resulting from 
inequality (fT3|). by comparing our estimation of (n), i.e. 

(n) (F) = ^ + - i , (20) 

with an independent measure of the mean number of photons. This comparison is discussed in the next section. 

We conclude this section by observing that if the state is not thermal, but we succeed in finding the value F 
saturating the uncertainty inequality to introduce an effective temperature T^ff. This effective temperature is again 
given by the right hand side of Eq. (jlSp and corresponds to a purity parameter tt (F). 

IV. EXPERIMENT 

Being the uncertainty relations related to optical tomograms through the tomographic function F (9) (see Eq. ((T2|) ) 
experimental data obtained by optical homodync detector, suitable for retrieving F {9), allow checking the uncertainty 
relation (fT3| . Moreover, F (9) allows to evaluate tt, the state purity (Eq. ([TT]) ). and in the case of thermal state, T, 
the field temperature (Eq. (fTS])). and (n), the mean photon number (Eq. ([20]) ). 

Thermal states are Gaussian so for these states also Eq. ([7]) is valid. Thus, in order to asses the reliability of the 
proposed method, it is possible to compare the estimations of tt via F (9) with the same quantity obtained by using 
Eq. ([7]). Moreover, a full reconstruction of the state via quantum tomography provides a further estimation of tt. 

To this end, F (9) has been retrieved for thermal continuous wave (CW) states, outing a sub-threshold non- 
degenerate optical parametric oscillator (OPO) 2^. In such a device non-linear fluorescence gives rise to a pair of 
down-converted entangled modes each in a thermal state [2l[. The experimental setup, illustrated in greater details 



elsewhere 



20 



2ll |. can be sketched into three distinct blocks: the state source, the below threshold OPO, the detector. 



a quantum homodync, and the acquisition board. 
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22, 1231 for 



The quantum homodyne detector, shows an overall quantum efficiency 77 = 0.88 ± 0.02 (see Refs. 
details). The system is set to obtain a 27r-wide linear scanning of the LO phase in an acquisition window. Since F (9) 
is retrieved by combining variances of data distributions calculated at different 9, we have decided to retrieve F (9) 
in [0,7r]. 

In order to use our homodync data to calculate F (9), we must be sure that the state is effectively a thermal one. 



First we prove that the state is Gaussian by performing some tests to assess the Gaussianity of data distribution [24 1 . 
In particular, we have used the kurtosis excess (or Fisher's index) and the Shapiro- Wilk indicator. The kurtosis is a 
measure of the " peakedness" of the probability distribution of a real- valued random variable while the Shapiro- Wilk 
one tests the null hypothesis that a distribution xi, ...Xn came from a normally distributed population. Then, a pattern 
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FIG. 3: F{e) vs LO phase 9 € [0,tv[. The mean values are: 0.65 ±0.01 (squares, tt = 0.48), and 0.329 ±0.007 (circles, = 0.61) 
where the errors are the standard deviations of the F {9) distributions in 9. The plotted range corresponds to 47 phase bins 
(see text for details). 

function tomographic analysis is used for assuring the thermal character of the state [l^ . [25^ and for reconstructing 
the state Wigner function. 

F (6) is calculated by analyzing the data distributions (each distribution containing 2100 data) at 47 different values 
of 9 (each 9 value then corresponds to a phase interval interval of 0.067 rad) and making use of Eq. 

A typical output is reported in Fig. [S] The traces are obtained by subtracting to each experimental value of F (9), 
Fshot corresponding to the average of F (9) over 6 for a vacuum state whose data are collected by obscuring the 
homodyne input. F (8) for a shot noise trace returns the of the instrument. Furthermore, in order to avoid any 
influence on the statistics of the data, the electronic noise is kept « 15 dB below the shot noise. The values of F (9), 
always positive as predicted by the uncertainty relation (jl3p . are with a good approximation independent of 9 as 
expected for a thermal state. 

We have analyzed 218 homodyne acquisitions of thermal states. This large number allows a statistical approach 
for analyzing the reliability of the proposed method. In particular the purity tt of the state has been evaluated: a) 
by reconstructing the state via homodyne tomography {ntom)', b) by using Eq. (|17p (ttf); c) by using the exact 
expression for a thermal state (see Eq. ([7])) (iTth)- In Fig. |4]we report the values of the differences Ap-tom (a) and 
Ath-tnm. (b) between fcp and ntom and iTth and iTtom normalized to their average. iTtom has been taken as a reference 
26|. By looking at these distributions it is seen that while Ath-tom is normally distributed (see inset in Fig. |4](b)) 
with an average of 0.022 (standard deviation 0.018) Ap-tom presents a sistemathic behavior in tt thus signalling the 
onset of systemathic error in the obtained determination. While nth comes from the exact expression of Eq. ([7]), TTi? 
is obtained by approximating (tt) (see Eq. (jlip ) so that the relative precision depends on the purity itself. It has to 
be noted that in any case the maximum relative error remains below 15%. 

Comparing np and nth evidenziate this effect even better (see Fig. [5]). 

The above discussed discrepancies, even not so small, confirms that the estimation of F{9)^ much more simpler 
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FIG. 4: Ap-tom (a) and Ath-to7n (b) vs. vftom- The histogram Ath-tom, given in the inset of (b), proves tliat Ath-tom is 
normally distributed. Ap-tom shows a well defined behaviour in ntom thus signalling the insorgence of systemathic error due 
to the used approximation. The estimation become more precise as the state approaches a pure one. The maximum error is 
arround 15%. 
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FIG. 5: Ap-th vs. TTtom- preseuts a systemathic behaviour in ntom due to the approximation of Eq. (fTTj) . 

than any full state reconstruction methods \l\ gives many reliable information on the liomodyned state. Moreover, 
for a thermal state it is possible to obtain more precise information by exploiting Eq. ([7]). 

By exploiting the thermal character of the states herein analysed it is possible to get information on T (Eq. ^TE\\ ) 
and (h) averaging the experimental trace for F (9) over 9. The values so obtained are given in the following table 
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(together with their confidence interval): 



mode 


F 


rpkB_ 


(h) 


1 


0.65 ±0.01 


3.8±0.1 


0.53 ±0.09 


2 


0.329 ±0.007 


2.8 ±0.1 


0.315 ±0.006 



(21) 



The analyzed mode has i' Ri 3 x lO^"' Hz so that the field temperature is of the order « lO'' K. 

V. CONCLUSIONS 

Uncertainty relations for mixed quantum states contain information that allows to easily estimate the purity 

once a connection between optical tomograms and these relations is set. In this paper it is shown that quadrature 
statistics of the field contains complete information on bound in the generalized uncertainty relations of mixed state 
and that this bound is intimately related to the state purity. This approach gives the possibility of a fast and reliable 
estimation of the purity parameter for quantum states that does not require sophysticate data analysis. As far as 
we know, all the methods for purity estimation rely on the reconstruction of the whole state density matrix while 
the approach presented in this paper allows to recover the purity simply calculating the tomographic function F (9) . 
The experimental results herein presented show that the homodyne detector can be used as the purity meter of the 
electromagnetic radiation in the quantum domain. The experimental estimations obtained via F (9) are compatible 
with more sophysticate estimation obtained via pattern function tomography so proving that our method can be used 
for real-time evaluations of some basic properties of the homodyne optical states. 

Further study of other photon states like squeezed states and multi-mode states will give the possibility to check 
other quantum phenomena as higher momenta uncertainty relations. 
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